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Abstract 

We study the effects on oscillations in the neutral-Kaon system from lepton and hadron 



\-{ ' self-energy loops in the VF-boson propagators. The corresponding W box diagrams are eval- 



uated by attaching the external quark lines to covariant vertex functions for the composite 
Kaons, integrating over all off-shell momenta. We find that the ratio of the imaginary and 
real parts of the amplitude is of the same size as the experimental value of the modulus of 
the CP-violating parameter e. 

The theory of weak decays, initially formulated by Fermi P1I2] in terms of a four-fermion 
vertex for the description of neutron (3 decay, is nowadays understood by a renormalizable jHlll] 
non-abelian gauge theory jSUnilZ], in which the weak interactions are mediated by the heavy gauge 
bosons. Here we shall concentrate on K^-K^ oscillations via double W exchanges. 

The phenomenon of a long-lived component in the two-pion decay mode of the neutral Kaon 
system, discovered by Christenson, Cronin, Fitch, and Turlay [H], is in the Standard Model (SM) 
parametrized by the complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) [HlCni matrix, 
which is related to the complex e parameter. This e measures the fraction of pion pairs observed 
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in the decay products of the long-hved neutral Kaons |lH ll2j. Many authors have studied this 
phenomenon, see e.g. Refs. P^TOITHIITMTIITMTIIIMTIO^^ iust to mention 

a few. But its dynamics remains to be discovered. 

The masses of the heavy gauge bosons, which can be estimated from Fermi's weak couphng 
constant, are experimentally determined with great precision jSU], ever since their first discovery 
in 1983 by the UA2 collaboration Furthermore, the branching ratios of the Ws and Z are 
well known for most of the decay modes that could play a role at Kaon energies. 

In scattering theory, unstable particles are described by Breit-Wigner amplitudes [S2], asso- 
ciated with complex-energy poles in the 5* matrix for the decay products Consequently, 
the physical, on-shell W boson, which besides being massive (mass M^), is also unstable (width 
Tw), should be characterized by the complex quantity Mw — ^iTw- 

In this paper, instead of describing the oscillations in the neutral Kaon system by a complex 
CKM phase and the exchange of massive W bosons with purely real mass Myy, we study whether 
a comparable result may be obtained by considering only real CKM matrix elements j33], but 
allowing the exchanged off-shell W bosons to develop a complex self-energy whenever this is 
kinematically demanded by four-momentum conservation and the known decay thresholds of the 
W. The crucial point here is that "decay is a profoundly irreversible process" [35j, which should 
also be taken into account in virtual exchanges, especially when hunting after a tiny effect. 
In order to have a reliable off-shell calculation, we shall sandwich the standard box diagrams 
between covariant composite-Kaon vertex functions that are parameter- free, after being tuned 
to the experimental Kaon size. Note that the Kaon's compositeness has recently been suggested 
to even give rise to CP violation [36 . 

In principle, one could calculate the desired amplitude by directly determining the two-pion 
branching ratio of the long-lived neutral Kaon ^7\. However, such a calculation would involve 
the unknows of QCD at low energies, which causes the result to depend on the specific modeling 
of strong interactions. Here, we determine the lowest-order processes that transform neutral 
Kaons into neutral anti-Kaons, schematically represented by the diagrams T and S of Figs. 
and respectively, and compare our result with the modulus of e. From experiment one 
derives \e\ = 0.0023 ± 0.0002. 

The strategy to obtain the modulus of e from a microscopic description stems from Gaillard 
and Lee [2^1 • It is based on a conjecture of Wolfenstein that all CP- violating processes 
stem from double-strangeness-changing [AS = 2) superweak interactions. The transition matrix 
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elements for neutral-Kaon oscillations are related to the modulus of e by 



(1) 



which equals half the ratio of the imaginary and real parts of the amplitude. 

The original calculus |39] was carried out by taking the external quarks free and the external 
momenta zero. However, we know that free quarks do not exist. Moreover, the neutral Kaons 
have long lifetimes — on a hadronic scale — thus justifying a bound-state approach with respect 
to their wave functions from strong interactions. Finally, the intermediate state (e.g. uu) may 
have a Kaon-like distribution w.r.t. strong interactions, but the involved u (or c, t) and u (or c, 
t) quarks can have all possible momenta, only governed by total four-momentum conservation. 
This picture is to be contrasted with the on-shell A/= 1/2 Ks — > tt+tt^ decay via an intermediate 
a resonance Hence, also for the intermediate state one must expect rearrangement effects in 
the quark-antiquark strong-interaction distributions. Therefore, we shall take the external quark 
momenta of the box diagram in agreement with the bound-state picture for neutral Kaons. As 
a consequence, our quarks will be dressed, and so massive. We take for the constituent quark 
masses the reasonable values m„ = iriu = md = 340 MeV and = 490 MeV [12] . In any case, 
we shall see later that the results hardly depend on these specific values. The initial and final 



quark-antiquark bound states couple, through the vertex functions, to the quark and antiquark 
participating in the ly-boson exchanges. The precise details of the vertex functions must be 
deduced from a microscopic model. We assume here that the vertex functions F and F depend 
on the relative quark-antiquark four-momenta p and p", respectively, and on the total center-of- 
mass (CM) momenta P = Pi + P2 resp. P" = pi" + P2" ■, i.e.. 



r 



r 




Figure 1: Diagram T: The t-channel box diagram for K - K oscillations. 



F = F(p,P) and F = F(p",P") 



(2) 



Total four-momentum conservation is expressed by 



P = P' = Pi'+P2' = P 



(3) 
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For the and vertex functions F resp. F we assume a Euclidean Gaussian distribution 
(exp — ), with a = 5.37 GeV~^, which in the nonrelativistic hmit corresponds to the 
experimental Kaon charge radius of 2.84 GeV~^ jHOlj and is, moreover, in good agreement with 
the predictions of the quark-level linear a model 1321 ^^^^ with the wave functions obtained 
in a unitarized meson model j44j . 

The relations between the particle momenta pi and p2 and the relative particle momenta p 
are given by the standard expressions 

pi,2 = ±p + "^"'^ P , (4) 
rrin + rus 

and similar relations for pi, p2 and p\ and pi", P2' and p" . The definitions Q have a nonrela- 
tivistic origin, but, as our calculation is fully covariant, the final result does not depend on this 
particular choice, which is nonetheless very convenient for numerical reasons. 

For the VF-propagators we write S^{qi, My/), where the exchange momenta qi and q2 are 
defined in Fig. [TJ according to qi = pi — pi and q2 = Pi" — Pi ■ The two-particle propagators 
Gqg{p, P) have the usual form 

and similarly for the other two. Now, we are dealing here with the propagators of quarks, which 
are fermions. Hence the two-particle propagators should be of the form 

^(^1,^2) = . . w . — -T ■ (6) 

[fi - mi le) \p2 - ^2 + ^e) 

This can be arranged by multiplying the two-particle propagators of Eq. (0) by 

(^1 mi) (^2 + rn2) . (7) 

However, for the crucial off-shell kinematical effects only the pole structure is of importance, 
which is the same for the fermion progators © and the boson propagators (0). Then, all 
the factors resulting from the Dirac algebra involving expressions like in Eq. ((7j), and also the 
different SM vertices, will cancel out when evaluating the ratio of the imaginary and real parts 
of the amplitude. Note that working with only scalar propagators allows us to use scalar vertex 
functions for the Kaons, and thus to model the diagrams in an essentially parameter-free way. 
The total amplitude for the process represented by diagram T is given by 

yl oc 1 dVdVc^VrG,,-5{^)G,^5£)G,jr. (8) 

The integrations are performed numerically, except for some angular integrations. 
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The pole structure of the two-particle propagators can best be studied in the CM frame (a/s 
represents the total invariant mass). If we then define 



n,s 



p +171 



2 

n,s 1 



(9) 



we obtain for the terms in the denominator of the two-particle propagator (0) the expressions 

2 



2 2 I • 



P2 - "^s + 



rrir. 



rUn + rus 



s+Po 



s-po 



uj„ + ie 



uj„ + ie . 



(10) 



+ rris 

We have poles in the complex pQ-plane for the two-particle propagator (0) at 

rrir. 



(1,2): po 
(3,4): Po 

graphically represented in Fig. 

(2) 

* 



m„ -|- m. 



■y/s ± (^uj'^ — iej and 



TTln + nis 



s ± iuj, — ie 



:iii 



(4) 
* 



3f?e (po) 



* 

(1) 



* 

(3) 



Figure 2: Singularity structure of formula ()1()|1 . The poles (*) are numbered according to their 
appearance in Eq. (fTTj) . 

Notice that < -|- rris in the bound-state regime. Therefore, the Wick rotation po —>■ ip^, 
and similarly for p' and p", can be freely carried out (see e.g. Ref. |13]), resulting in a purely real 
amplitude, at least when neglecting the W self-energy effects. 



5 



For the full determination of (k^ K^^, we also include the s-channel double-Ty-exchange 
diagram depicted in Fig. El We obtain contributions for diagram 5" which in magnitude are com- 
parable to those for diagram T. Moreover, diagram S requires less numerical effort than diagram 
T, since in the former case all angular integrations can be performed analjd;ically. Consequently, 
checking details on how our results depend on a and the quark masses we only do for diagram 
S. 




Figure 3: Diagram 5": The s-channel box diagram for - oscillations. 

Now, let us finally come to the main issue of this paper. In order to treat the exchanged 
W resonances in a more realistic way than what is normally done, we first include self-energy 
bubbles jiB] corresponding to the well-known [30] leptonic decay modes of the W , i.e., ei/g, /xz/^, 
and Tfr, as depicted in Fig. HI 

+ --<z>-- + - 

Figure 4: The full W propagator dressed with lepton loops and possibly with hadron loops 
Hence, we substitute in the W propagators by M^r + + + Z^, according to 

= -- 1-^ , 12 

2 \Pi\ \ r/ V PMink./ 

where F^, me {i = e, /x, r) and are taken from Ref. [SD]- Note that in the 9 function 
one must take the Minkowskian = pi — p"^, with po real, too. For obvious reasons, we have 
only kept the imaginary parts of the contributions to the mass self-energy (Fig. H}, the real parts 
being absorbed in the physical W central mass. At this point, we should also stress that adding 
a negative imaginary part (fT^ to the W mass does not spoil our straightforward Wick rotations, 
as it only amounts to substituting an infinitesimal e by a finite |Z(p^)|. In other words, because 
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of causality the sign of the imaginary part of the mass is always negative, both for stable and 
unstable (anti)particles jlTj. 

In Tabled we show how much each of the terms Zi contributes to the ratio in Eq. (Q). 



mode 



CS (DK) 

cX (DB) 



fraction 130 



(10.72±0.16)% 
(10.57±0.22)% 
(10.74±0.27)% 

(3111?)% 
(33.6±2.7)% 



1 






3?e( 







0.00096 
0.00094 
0.00060 
0.00153 
0.00078 



Table 1: Contributions to the ratio of the imaginary part over the real part of the K'^-K^ 
amplitude from the most important W decay modes (see also text). 



For the hadronic modes we make the following considerations. The only known mode con- 
taining a light hadron is tcj, having a decay rate of less than 0.1% [30], hence negligible. The 
next level of known modes contains at least one c quark. The threshold for the DK loops at 
2.37 GeV is some 600 MeV higher than for tUt, whereas the branching ratio is about the triple 
of the latter, which yields a contribution about two and a half times larger than from (see 
Tabled). 

The fraction of cX decay modes of the 14^ is a bit larger than three times the fraction of 
tUt-. However, here we do not have much information on X. If we assume X = b, we get a 
contribution from DB loops which is smaller than from DK, but still considerable. 

It is clear from these considerations that we cannot include the hadronic loops to a great 
accuracy, because of the experimental uncertainties. Nevertheless, we may conclude from our 
results that their total contribution will not be very different from the sum of the two individual 
cases included by us, which is a kind of upper bound. Nevertheless, we shall take a conservative 
choice for the total final theoretical error in our result. Summing the contributions of Table Q 
(we actually put all loops together when determining the total (K^ K^) amplitude), we obtain 



(13) 









23fte ( 







0024 

U.UUZ^ -0.0004 



for the total contribution of W self-energy loops, which is of the same magnitude as \e\ 

In the errors quoted in Eq. (fT^ . we have included, besides small numerical errors, the effects 

of reasonable variations in the parameters a, rrin, and m^, as well as the uncertainty in the 
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cX mode. Also the contributions for intermediate c and t quarks have been studied, which do 
not significantly alter the result (fT^ . though adding to the total (^K^ K^^ amplitude. The 
coincidence of the effect computed by us and \e\ is very striking, the physical significance of 
which is open to debate. 
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